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ABSTRACT 



Aims. In the dense stellar environment of the globular clusters, compact binaries are produced dynamically. Therefore the fraction of 
type la supernovae that explode in globular clusters is expected to be higher than the fraction of mass residing in these. 
Methods. We have searched for globular clusters at the positions of observed type la supernovae. We used archival HST images 
and literature data, covering the positions either before the supernovae exploded, or long enough after that the supernovae have faded 
below the luminosities of globular clusters. 

Results. We did not find evidence for globular clusters at any of the supernova positions. For 18 type la supernovae, the observations 
are sensitive enough that any globular cluster would have been detected, and for further 17 type la supernovae, the brighter globular 
clusters would have been detected. Correcting for incompleteness, we derive a 90% upper limit of 0.09 on the fraction of type la 
supernovae that explode in globular clusters for the full sample and 0.22 for the sample of supernovae in late-type galaxies. This 
allows us to limit enhancements per unit stellar mass for a coeval population rj^o ^ 50 (100) with 90% (99%) confidence. We find that 
by observing the positions of a sample of less than 100 type la supernovae in the outer parts of early-type galaxies, it will be possible 
to probe the currently favoured range of t/^o ~ 1 - 10. 

Key words. Supernovae: general ~ Galaxies: star clusters: general 



1. introduction 

Type la supernovae (SNIa) are believed to be thermonuc l ear ex - 
plosions of white dwarfs (e.g. iHillebrandt & Niemeveii l2000l) . 
The two main scenarios for bringing the white dwarfs above 
the critical explosion mass Mc (similar but not equal to the 
Chandrasekhar mass) are (1) the single-degenerate (SD) sce- 
nario, in which a white dwarf accretes from a non-de generate 
companion star (IWhelan & Ibenill973UNomotolll982h . and (2) 
the double degenerate (DD) scenario, in which the merger of 
two white dwarfs with a total mass exceeding Mc causes the ex- 
plosion (Iben & Tutukov, 1984; Webbink, 1984). Both scenarios 
requires binaries with white dwarfs in orbits close enough that 
mass transfer occurs at some point during the evolution of the 
binary. 

Globular clusters (GCs) are known to harbour large popula- 
tions of compact binaries. Most well known are the bright low- 
mass X-ray binaries (LMXBs) with a specific density per unit 
mass 100 times higher than in the field (Clark, 1975). Also the 
popul ation of millisecond pulsars (Lvne et al.. .1987; Lorimerl 
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l2005h and blue stragglers (ISandagel Il953h have been found to 



be strongly enhanced. Particularly interesting for type la super- 
novae are the relatively recent observational evidence that also 
the populations of white dwarf systems are enhanced (e.g. Hei ; 
Heinke et al.1 120051: iDieball et all l2007t iMaccarone & Knigge . 



2007tlKnigge"eralTl2008tlHenzeetal.Ll2009ir The enhancement 
of tight binaries in globular clusters is a consequence of the very 
high stellar densities found there (up to IO^Mq pc"^). With such 
high densities, dynamical encounters which create or modify bi- 
naries are frequent, and the higher mass of the binaries makes 
them sink to the center, where the encounter rates are highest. 

It is therefore reasonable to expect that both SD and DD 
SNIa progenitors will be enhanced in globular clusters. The 



few theoretical studies that have investigated this have indeed 
found such enhancem ent in their models (Ivanovaet al., 200^ 
IShara&Hurievl 120061) . However, the magnitude of such an en- 
hancement is highly uncertain. Aside from the fact that it is not 
known which systems actually do lead to supernova explosions 
and how these systems evolve outside globular clusters, there is 
also a very large uncertainty due to the difficulties of modelling 
and observing exotic binary populations in globular clusters. 

The most massive globular clusters, with the highest colli- 
sion rates are bright and can therefore be seen to large distances 
(~50-100 Mpc). It is therefore feasible to put limits on the frac- 
tion of type la supernovae that explode in globular clusters by 
identifying these in deep pre- or post- explosion images of the 
supernova positions dPfahl et all l2009h . Despite this prediction, 
no such observational survey has been carried out. 

In this paper we perform a survey using archival observations 
and literature data to place observational constraints on the frac- 
tion of SNIae in globular clusters. We first discuss the theoretical 
estimates in Sect. |2] then we discuss our methodology in Sect. 
[3] In Sect. |4] we describe and analyse the literature and archival 
data, and in Sect. |5] we discuss the results and perspectives for 
future surveys. 



2. Theoretical expectations 

It is clear from both observations and theory that compact bina- 
ries can be formed efficiently through dynamical interactions in 
the dense stellar environment of globular clusters. However, the 
picture of the formation and evolution of these compact binaries 
remains vague, due to a large number of unconstrained physical 
processes involved, and the computational expense of globular 
cluster models including binary evolution. Modelling the forma- 
tion of SNIa progenitors is even more problematic, given their 
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unknown nature. It is therefore not surprising that only few stud- 
ies have attempted this. 

The enhancement of compact binaries in a globular cluster 
is expected to be proportional to the rate of interactions F (e.g . 
iHut & VerbunR[T983[|Poolev et all I200I iPoolev & Huil2006l) . 
While r can be estimated for nearby globular clusters, based on 
their structural pa rameters, the estimates are not very reliable 
dMaccarone et al.L |201 1). and it is not possible to measure the 
parameters accurately outside the Galaxy. For this reason we use 
the commonly used average enhancement factor for a galactic 
population of globular clusters 



2: Mccj X T]cc,i 

N 

Z Mccj 



(1) 



Single-degenerate p rogenitors were considered in the study 
of llvanova et al.l (l2006l) . who found an enhancement factor (per 
unit stellar mass) of 77= 1-7 compared to a field population with 
solar metallicity. They find no single-degenerate SNIa in a pop- 
ulation of stars with the same properti es as the globul a r clus- 
ters but with interactions turned off". IShara & HurlevI (l2006l) 
finds a small enhancement <2 in the production rate of accret- 
ing white dwarfs, but find that in their globular cluster models, 
the accreting white dwarfs are heavier, and so more likely to 
be SNIa progenitors. Other stu dies have predicted higher num- 
bers of accreting white dwarfs dPi Stefano & RappaportL Il994t 
iDavies & BenzlTr995b . correspond ing to ? 7=5, but did not con- 
sider SNIa progenitors specificallv. IShara & Hurley (2006) find 
no enhancement of double white dwarf binaries, but this study 
does not discuss SNIa progenitors specifically, llvanova et alj 
(l2006l) find an enhancement of merging double white dwarfs 
above the Chandrasekhar mass of 77 =3-13 compared to a field 
population with the same metallicity and age as the popula- 
tion in the globular clusters, but no significant enhancement 
(77 - 0.3 - 1.2) when compared to field population with solar 
metallicity. 

For the SD scenario, an alternative approach is to compare 
to observations of similar systems. This is not possible for the 
DD scenario, due to a lack of observational systems to com- 
pare to. With their accreting massive white dwarfs, SD SNIa 
progenitors are somewhat similar to cataclysmic variables (CVs) 
and LMXBs in their formation and evolution. The bright LMXB 
population in globular clusters is quite well studied, as they can 
be seen to large distances with Chandra, and they have been 
found to be over- abundant by a factor of -100 ( C laiS, Jl975t 
LSarazin et al., 2003; Jordan et al., 2007; Voss et al., 20091^. The 
CV population is much less understood, as they are much harder 
to identify. Only a small sample has been found in recent years, 
consistent with an over-production by a factor of ~few (e .g. 
iPoolev & Hut. 2006"; Dieball et al., 2007; Knigge et al.', 2008), 
but completeness is a serious issue, and higher enhancement fac- 
tors are therefore not ruled out. Obser vations of novae i n M3 1 
suggests en enhancement factor of ~ 10 dHenze et alil2009l) . The 
SD SNIa progenitors have white dwarfs with masses near Mc, 
more similar to the masses of neutron stars than to those of most 
CVs. They therefore sink the the center more easily and experi- 
ence more dynamical encounters than more typical white dwarf 
systems, which leads to a higher expected enhancement. 

From the discussion above we conclude that 77 is most likely 
greater than one and lower than 10. However, it is clear that the 
results are very poorly constrained, and models with 77 outside 



this range can not be discarded. It is therefore important to con- 
strain this fraction from observations. 

Besides 77, the fraction of SNIae that explodes in globular 
clusters also depends on the fraction of stellar mass that resides 
in the globular clusters, Fm,cc - Mcc/Mf, where Mcc is the 
total mass of the population of globular clusters and Mf is the 
total mass of all the other stars. For a sample of SNIae, the 
expected number that explodes in globular clusters is then: 



Ncc ^Nx Fm.cc X 77. 



(2) 



The mass fraction Fm,gc varies strongly between galaxies. The 
Milky Way has a low abundance o f globular clu sters, with 
Fm.cc ~ 0.1% (e.g. the catalogue of iHarrisl [19961) . Therefore 
only ~ 10% of the Milky Way LMXBs (see e.g. iLiu et al.L 
120071) . and based on the considerations above, probably ~ 1% 
of the CVs, are in globular clusters. This is in stark contrast 
to many elliptical galaxies with rich globular cluster s ystems, 
where in some cases Fm,gc can be larger than 1% (e.g. Harrisl 
2009)). and the majorit y of bright LMXBs are in globular clusters 
( Angelini et all 1200 lb . Correspondingly, if the SNIa enhance- 
ment is a factor of 10, about 10% of the SNIae in these galaxies 
must be formed in globular clusters. Averaging over the popu- 
lation of nearby galaxies yields a fraction of SNIae in globular 
clusters of ~ 1 - 3 ■ 10"^?7 (Pfahl et al., 2009), most likely a few 
per cent. 

A complicating factor that was not taken into account in pre- 
vious studies is the fact that the rate of type la supernovae RsNia 
decreases with age for a coeval population of stars. In general 
the distribution of globular cluster ages is different from that of 
the field stars. Therefore 
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RsNIa(tcc) 



(3) 



where tec is the age of the globular clusters and f/ is the age 
of the field stars. For early-type galaxies tec ~ tf, whereas the 
bulk of the field population of late-type galaxies tends to be 
significantly younger than the globular cluster population. 77^0 
is the globular cluster enhancement per unit stellar mass, com- 
pared to a field population of the same age. The exact shape of 
the delay-time distribution (DTD, the SNIa rate as a function 
of time for a coeval population of stars) is not known, but it 
has been shown to decrease strongly by a factor of > 10 from 
youn g environrn ents <1 Gyr to older environments ~ 10 Gyr 
(e.g. lMaoz et"!!]. 12010). Therefore -50-85% of all type la su- 
pernovae explode within the first Gyr after star formation (e.g. 
iMaoz et aTL I2OIOI) . However, the current star-formation density 
is much lower than at redshifts 1, and comparing the local rate 
of star formation |Hanish et al. (2006) to the stellar mass density 
Salucci & Persic ( 1999): lCole et al.1 (12001) leads to a fraction of 
stars in the local universe that were formed less than 1 Gyr ago 
of ~ 2%. Despite the high fraction of prompt SNIae for a coeval 
population of stars, the local population will therefore be dom- 
inated by the tardy component. Combining the fractions found 
by Sullivan et al. (2006) with the local rate o f star formation 
Hanish et al.. (.2006.) and stellar mass density iSalucci & Persid 
(ll999h ; ICole et al.1 ( 1200 lb . only < 20 per cent of the local SNIae 
are expected to belong to the prompt component. 

will 



For most local galaxies, the ratio Rgc/f = 



therefore be higher than what could be expected from the DTD. 
It is typically in the range ~ 0.1 - 1. For early-type galaxies, 
the typical age of the field population is similar to the age of 
the globular clusters, and therefore Rcc/f ~ 1- Late-type galax- 
ies can have significant populations of young stars for which the 
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SNIa rate is more than a magnitude higher than for the old pop- 
ulation of stars in their globular clusters. However, the vast ma- 
jority of late-type galaxies also have older stellar components, 
and Rcc/F will therefore almost always be higher than 0.1. 



3. Exclusion of cluster origins for a sample of type 
la supernovae 

The association of type la supernovae with globular clusters re- 
lies on the spatial coincidence. If a type la supernova is found to 
have exploded at the same position as a globular cluster, it may 
have exploded inside the globular cluster If not, then a glob- 
ular cluster origin is definitely excluded. However, the distri- 
bution of GC luminosities is wide, and many observations are 
only sensit i ve to the bright end of the GC luminosity function. 
iPfahl et al.1 (|2009|) showed that for 77= 1 approximately 1 % of all 
SNIa should explode in globular clusters, corresponding to 1 su- 
pernova in a globular cluster within 100 Mpc every year, and that 
a dedicated HST program would be able to find the connection, 
if each SNIa was observed a few years after the explosion. With 
their assumptions such a program would require deep HST ob- 
servations of ~ 100 fields every year They did not account for the 
fact that some SNIa suffer from significant extinction and that 
observations inside galaxies where there is a high background 
are much less sensitive than in the field. These effects are diffi- 
cult to model for the full sample of type la supernovae, and it is 
therefore unclear if the proposed observational program would 
be succesful despite the high costs. 

We therefore use the currently available HST data to observe 
or put limits on the fraction of SNIae in globular clusters, and to 
make it possible to estimate the sensitivity that can be achieved 
with a dedicated observational program. 

3.1. Globular cluster completeness 

To be able to exclude a GC origin, it is necessary to know the ob- 
servable properties of the GCs. Being old stellar systems, there 
are only relatively small variations in the mass-to-light ratios 
of different clusters, with the main difference being related to 
the globular cluster metallicity. However, the mass distribution 
of globular clusters is wide, with several orders of magnitude 
difference between the brightest and the faintest clusters (e.g. 
[Harris, 1991; Jordan etal., 2007). In the more distant galaxies 
or in shallow observations it will not be possible to observe the 
faint clusters and it is necessary to calculate the incompleteness 
caused by this. 

The most straight-forward way to do this is to use the 
observed luminosity function of globular clusters dPfahl et al.L 
I2OO9I) . However, while the definition of rj as the enhancement 
per unit mass is a useful measure, it is misleading in terms of the 
physical interpretation, as the number of compact binaries does 
not scale linearly with the mass of the globular cluster Instead 
it is related to the stellar encounter rate F inside the clusters. 
The distribution of structural parameters (and hence F) of clus- 
ters is not well known. Since also the exact processes of compact 
binary formation are poorly constrained, it is therefore not pos- 
sible to make theoretical estimates of the relation between GC 
mass and the probability of hosting SNIa progenitors. 

We instead attempt to do this on an empirical basis. The only 
compact binaries that have been surveyed in large samples of 
GCs are LMXBs. We us e the results of Peacock et al. (2010) 
and ISivakoff etaP JIoOTI) to estimate the mass distribution of 
GCs that contribute to the population of compact binaries, and 



Table 1. Empirical completeness levels (absolute magnitudes) 
of globular clusters with LMXBs in the K- and z-band, and the 
corresponding globular cluster masses (in 10^ M©) from simple 
stellar population models. 



Completeness 


CLIOO 


CL75 


CL50 


CL25 


K-band 


-9.5 


-11 


-12 


-12.5 


z-band 


-7.5 


-9.0 


-10.0 


-10.5 


Mz=0.012 


1.7 


6.7 


17 


27 


M z=0.0012 


2.5 


9.9 


25 


40 



thereby to estimate 4 completeness levels, CL25, CL50, CL75 
and CLIOO, meaning the GC masses/luminosities above which 
25%, 50%, 75% and 100% of the compact binaries are expected 
to reside. We strongly caution that this is based on observations 
of LMXBs alone, and that it is very possible that SNIa progeni- 
tors could have a different mass-dependence. However, currently 
there are no models that predict a different behaviour, and the 
populations of faint X-ray sources in Galactic GCs do seem to 
follow the distribution of bright LMXBs. 

To estima te the com pleteness masses, we use the K- 
band data of iPeacock et al.l (120101) and the z-band data of 



ISivakoff et al ] J2007h . We give the completeness values for the 
K-band observations of M3 1 and z-band observations of Virgo 
in table [T] To compare these results and extrapolate them to 
other wavelengths we use the int egrated simp le stell ar popula- 
tion magnitudes of iGirardi et alJ I2OOOI) : iMarigo et alJ (i2008h . 
The values are in good agreement, assuming a 12 Gyr stellar 
population with a Chabrier initial mass function and a metal- 
licity of 0.012. We also provide estimates of the corresponding 
(initial) globular cluster masses, using the K-band magnitudes 
and two different metallicities. These are then used to find the 
colours in all the different bands used in the following analy- 
sis. As the metal-rich GCs are redder, they are fainter than the 
metal-poor GCs in the used bands, for a given K-band luminos- 
ity. We therefore use the calculations for z=0.012 to determine 
the confidence limits in all the bands, noting that in this way the 
confidence limits will be underestimated, decreasing our sensi- 
tivity somewhat. Younger clusters would be brighter for a given 
stellar mass, similarly leading to an underestimation of the con- 
fidence limits. The magnitudes for CLIOO are given in table |2] 
From table[T]it can be seen that subtracting 1.5, 2.5 and 3.0 from 
these magnitudes yields CL75, CL50 and CL25, respectively. 

4. Results 

The observations of globular clusters at the positions of type la 
supernovae can be performed before the explosions, but also af- 
ter, since the light of the globular clusters is provided by stars, 
which are not affected by the explosion. As the luminosity of the 
SNIae are very high just after the explosion, it is necessary to 
wait until they have become fainter than a globular cluster This 
typically takes less than two years, although SNIae with light- 
echoes might be bright enough to hide globular clusters for a 
longer time. 



4.1. Literature survey 

To facilitate comparisons and thereby the use of type la super- 
nova as standard candles in cosmology, their lightcurves are ob- 
served using BVRI photometry. While most studies are limited 
to less than ~ 100 days, a number of SNIae have published late- 
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Table 2. Assumed absolute magnitudes for CLIOO. The HST magnitudes are for the VEGA photometric system. The differences 
between the magnitudes with different instruments is too small to be of importance for our study. 



Band 


B 


V 


R 


1 F555W 


F569W 


F675W 


F814W 


F850LP 


F160W 


F190N 


Mag 


-5.5 


-6.5 


-7.1 


-7.7 -6.5 


-6.6 


-7.2 


-7.6 


-7.9 


-9.3 


-9.5 



Table 3. Faintest observed BVRI magnitudes or upper limits from late-time observations of type la supernovae in the literature. CL 
gives the completeness of ruling out globular clusters at the position, see text. 



Supernova 


B 


V 


R 


1 


Galaxy 


Gal type 


Distance 


Av 


References 


Type CL 


SN1990O 


22.747 


22.886 


23.629 


22.368 


MCG+03-44-03 


Sba 


150Mpc' 


>0.29"' 


8 





SN1990N 




23.198 






NGC 4639 


Sb 


27.0Mpc' 
18.5 Mpc- 


0.221' 


9 


25 


SN1991bg 




25.0 






NGC 4374 


E 


0.096' 


10 


P(sublum) 100 


SN1992A 


>26.5 


>26.5 






NGC 1380 


SO 


21.2Mpc- 


0.014' 


10 


100 


SN1992bc 


22.716 


22.190 


23.172 


22.935 


ESQ 300-09 


s 


lOlMpc' 


0.012' 


12 





SN1993L 




23.0 




21.5 


IC 5270 


So 


23.7Mpc^ 


> 0.23"^ 


10 





SN1996X 


23.68 


23.89 


21.63 


20.78 


NGC5061 


E 


25.5Mpc^ 


0.031' 


11 


25 


SN1997cn 




>23.198 






NGC5490 


E 


78.3Mpc' 


0.0'^ 


12 





SN2000ce 


24.12 


23.84 


23.77 


22.84 


UGC4195 


SBc 


86.7Mpc' 


1.67' 


13 





SN2001C 


22.64 


22.89 


23.57 


23.23 


LEDA 19975 


Sb 


47.6Mpc-' 


0.403" 


13 


25 


SN2001V 


22.60 


21.96 


23.09 


21.84 


NGC 3987 


Sb 


68.5Mpc' 


0.171' 


13 





SN2001bg 


21.99 


22.10 


22.61 


21.61 


NGC 2608 


SBb 


36.3 Mpc- 


0.868" 


13 





SN2001dp 


20.40 


21.49 


21.38 


19.62 


NGC 3953 


SBb 


n.lMpc^ 


0.09"' 


13 


25 


SN2003du 


22.771 


22.827 


23.010 


22.121 


UGC 9391 


SB 


44.5Mpc' 


0.032' 


14 





SN2006gz 


>24.4 


>24.2 


25.5 




IC 1277 


So 


103Mpc-' 


0.753" 


15 


P (overlum) 


Sandase et al. 


(2010),: 


■ Tiillv et al. 1 200^, ^Average redshift-independent distance from NED, 


redshift-based distance from NED, ' Hicken et al. 1 2009), ^ Cai 



T\Wk . iKnsciunas et al.gOOU . lHamuv et al.a996i) . iLrra et aljfT99l . " jMiTne et al. 1 1999i).' 1Saivo et aUtlOOil) , ' ITuratto et al] fT99l , ' iLair et al., a006.) , '^ .Stanistiev et ■dUamxk . 

' IMaeda et alH200l , " ISchlegel et aj] (1993) . 



time data. We have surveyed the litterature to compile a sample 
of 19 such SNIa. None of the published results show signs of a 
constant component expected from a host cluster. We list these 
in table[3] The table also lists the faintest observed magnitudes in 
the BVRI bands, as well as the distance to the host galaxy and an 
estimate of the extinction towards the supernovae. We use these 
to find the intrinsic absolute magnitude of the faintest observa- 
tion in each band and compare these with the values given in 
table [1] and |2] to find the confidence with which we can exclude 
globular clusters on the basis of these light-curves. In table|3]we 
have left out four SNIae (SN1991T, SN2000E, SN2000cx and 
SN2001el) for which useful late-time lightcurves exist, as more 
constraining limits are found in the analysis below. 

In table |4] we have compiled a second more heterogenous 
literature sample. This consists of two old SNIae that were ob- 
served before the definition of the BVRI photometry (SN1937C 
and SN1972E), one where the magnitude was reported using 
HST with the VEGA magnitude system (SN2000cx) and 5 
where upper limits were published based on pre-supernova HST 
images. Similar to table |3] we derived completeness limits for 
these supernovae. 

4.2. HST observations 

In addition to the literature survey, we have searched for archival 
HST data at the positions of all known SNIa within 100 Mpc. 
We analysed the data to either find observations where the su- 
pernova is faint enough to exclude a GC origin, or where the su- 
pernova was not observed and an interesting upper limit could be 
inferred. One problem is the positional accuracy of the sources, 
as the fields of external galaxies can be crowded. Therefore po- 
sitions of sub-arcsec precision are needed. Unfortunately many 
supernova positions are relatively poorly determined, as the only 
published coordinates are from the original detections with small 
telescopes, when the supernovae are very bright, and the coordi- 
nates are provided without error estimates. In the case where 



several groups have published coordinates, the distance between 
the positions can be large. Unfortunately most supernova obser- 
vations are not publicly available, making it impossible to verify 
the positions. 

We therefore only report results from supernovae, where we 
are confident that the positions are well-determined. In a few 
cases, the limits were obtained from images, where the super- 
nova is still seen at a luminosity below that of globular clus- 
ters. In some cases, the positions could be found from other HST 
images taken when the supernova was visible, which were then 
matched to the image from which the limit was determined. In 
the cases where the positions could no t be determined from HST 
images, we used the positions listed in lHicken et al.l (l2009l) . who 
performed careful astrometry of a large sample o f supernovae. 
For ea ch of the supernovae where we used the ' Hicken et al.l 
(l2009h positions we compared to HLA HST images, where the 
HST astrometry was already corrected using either 2MASS, 
GSC or SDSS as reference. We also used these catalogues to 
confirm the astrometry, as the HLA astrometry is done automat- 
ically. We furthermore searched for ground-based observations 
of the supernovae, which allowed us to find the positions of 7 
additional SNIae. 

To find the upper limit in an individual observation, we used 
the background counts in a number of circular regions near 
the supernova position (excluding point sources) to estimate the 
count rate of a 3-sigma fluctuation C30-. The upper limit Cul 
on the supernova counts were then calculated by subtracting the 
median of the background fields Cul = C30- - C^. The ra- 
dius of the circular regions were chosen in the range 0.15-0.5 
arcsec, depending on the local density of point sources and the 
gradient of the host galaxy light. From this we calculated the 
upper limit in the VEGA magnitude system using the standard 
count-rate to magnitude conversions from the HST data hand- 
book, including PSF corrections for the aperture size. We teste d 
our method against the method used in iNelemans et al.l (l2008l) : 
IVoss & NelemansI (l2008i) . where fake sources were put in the 
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Table 4. Additional published optical limits on SNIa in GCs, with indication of whether the limit was from images obtained before 
or after the supernova explosion. HST magnitudes are in the VEGA photometric system. CL gives the completeness of ruling out 
globular clusters at the position, see text. 



Sup Grn o V3 


When 


Opticfll result 


Gtiltixy 


Gtil type 


Distance 


A,, 


References 


Oi 'I 1 J 3 I 


After 


\/T'j nni tiifli=> 90 


IC 4182 


Sdm dwiirf 


4 IMpc" 


> 0.04"' 


21 


SN1972E 


After 


Magnitude 17 


NGC 5253 


10 


3.4Mpc- 
24.0Mpc- 


> 0.17" 


22 


SN2000CX 


After 


F555W=25.2 


NGC 524 


so 


0.248^' 


23 


SN2003cg 


before 


F814W=22.9 


NGC 3169 


Sb 


21.6 Mpc' 


> 0.10"5 


32 


SN2003gs 


before 


F555W=25.1 


NGC 936 


SBO 


23.0 Mpc- 


0.205'* 


32 


SN2003hv 


before 


F814W=25.7 


NGC 1201 


SO 


20.2 Mpc^ 


0.05" 


32 


SN2004W 


before 


F850LP=26.0 


NGC 4649 


SO 


16.5Mpc' 


2.4™ 


32 


SN2005df 


before 


F606W=26.8 


NGC 1559 


SBc 


15.4 Mpc-* 


< 0.40^* 


32 


SN2006dd 


before 


B=26.0,V=26.1,I=25.4 


NGC 1316 


Sa 


17.8Mpc-° 


0.25^* 


24 


SN2006mr 


before 


B=25.8,V=26,I=25.4 


NGC 1316 


Sa 


17.8Mpc-° 
22.0Mpc- 


0.25" 


24 


SN2007sr 


before 


F8I4W=25, F555W=26.5 


NGC 4038 


Sc 


0.558" 


25 


SN20{)7on 


before 


F475W=27 


NGC 1404 


E 


20.2Mpc2 


> 0.03'" 


26 


SN20{)8ge 


before 


F606W=24.3 


NGC 1527 


SO 


18.0Mpc^ 


0.040-' 


27 


SN2011fe 


before 


F814W=26.35 


M81 


Sb 


6.4 Mpc^* 


> 0.04'" 


29 



Type 



CL 



P(91T-like) 



P(sublum, fast-decl) 



P(02cx-like) 






75 
75 
75 
100 
100 
100 
100 
100 
100 
100 
75 
100 



References; 1-16 see table[3] ' 'Milne et al. i201Q.).' ^iKn"sciunas et aljt2009l) .' lLeioudas et alj t2009h .^ 1Stritzin ger et al.l>20T^.^'Baade et aljfT953). ^-"Kirsh ner &"Okg fT97l . 
- ISoUerman etaj] 120041) . ^fMaoz & Mannucci 1 200l . ^ INelemans eTS] 12003) . ^ 1Voss & Nelemanj l200l . ^ ' Foley et al. 1 20 1 0), Shappee & Stanek ( 20 1 1 ), - lLi et all t20Tll) . 
^ 1Elias-Rosa et alJI20Q6l) .^ lLi et aljt2001i) . ^"S.J. Smaitt, private communication, see |http : //www, lorentzcenter .nl/lc/web/201SI/391/presentations7Smartt -pptj 



images and the upper limit was based on the detection of these, 
and we found good agreement. 

The results are listed in table|5] for the supernovae where the 
upper limits provide constraints on the globular cluster connec- 
tion. 

4.2.1. Notes on individual sources. 

Some of the analysis of individual sources needs more explana- 
tion. SN1998bu and SN1991T have light echoes that are seen 
in all the HST observations. For these sources we use the last 
F814W observations, and we measure their flux within an aper- 
ture of 0.2 arcsec radius from the central source position. For 
both SNIae, the flux from this region is low enough that the ex- 
istence of globular clusters at these positions can be excluded. 

For many of the SNIae in table |5] the position can not be 
found from HST data. Four positions (S N2002bo , SN20 02fk, 
SN2005el and SN2005hk) are from Hic ken et al.l (l2009h . and 
further 7 were found from ground-based data that we matched 
to the HST images. Four SN1994ae, SN2001el, SN2005df 
and SN2007af were identified in images taken with the NTT 
(SN1994ae) and VLT (SN2001el, SN2005df and SN2007af) 
obtained from the ESO archive. Three more SNIae positions 
were found from ground-based images provided by Weidong 
Li. SN1999gh and SN1999gd were observed with KAIT and 
SN1998aq with a 1. 2m CfA telescope. 

4.3. Derived upper limits 

The results above constitute the first observational survey of the 
connection between globular clusters and type la supernovae. 
From this it has been possible to exclude a connection fully or 
partially for 35 SNIae, and no SNIa with possible globular coun- 
terparts has been found. Taking into account that faint globular 
cluster counterparts would still be possible for some of the SNIae 
in our sample, this allows us to derive Pcl - 90% and 99% up- 
per limits on the fraction of SNIae in globular clusters Fm. If a 
fraction FsNia of SNIae explode in globular clusters, the prob- 
ability of not finding a globular cluster at the position of SNIa 
number / is 1 - Fs Nia x CL,, where CL, is the globular cluster 
completeness level (which for our sample can take the values 
0.25, 0.5, 0.75 and 1.0, see section [TTT l. The probability of find- 
ing no globular clusters for the entire sample of SNIae is the 



product of finding no globular clusters for each SNIae and Fj/l, 
the upper limit on FsNia, can therefore be found by solving 

N 

Y](l-FuLxCLi)^l-PcL. (4) 

!=I 

Table |6]lists the limits for the full sample, and for two sub- 
samples. The first subsample consists of all the normal SNIae, 
which follow the same standard lightcurve evolution. This sam- 
ple is chosen to eliminate the ones with peculiar lightcurves, as 
they are more Ukely to have been studied in detail and therefore 
can bias our sample by being over-represented. Furthermore, 
their origin might diff'er from the normal SNIae. We note that the 
peculiar label is subjective and some of them could belong to the 
normal population, making the normal sample overly conserva- 
tive. Our second subsample consists of all the SNIae in early- 
type (elliptical and SO) galaxies. These have old stellar popula- 
tions that are similar to the populations found in globular clus- 
ters, and are therefore particularly useful for comparison. 

The 90% upper limit from the full sample is very close to 
our expectation that a few per cent of type la supernovae could 
explode in globular clusters, and it is therefore useful for con- 
straining theoretical models. From equations |2] and |3] the upper 
limits on Fm corresponds to upper limits on 77: 



Fm.gcX < Rcc/F > 

For the full sample, the average stellar population is significantly 
younger than in the globular clusters, and we therefore assume 
< Rcc/F >~ 0.2-0.4. As the fraction ofglobular clusters is small 
in the late-type galaxy part of the sample, Fm,cc ~ 0.1 - 0.3%. 
With these assumption we get a result of rfco < 75 - 450. For the 
sample of early-type galaxies, the ages of the stars are similar to 
the ages of the globular clusters and therefore < Rgc/f >~ 1- 
Furthermore, the globular cluster mass fraction is higher, and 
we assume an average value of 0.3 - 0.6%. This leads to a more 
constraining limit of rjco < 37-73. 

5. Discussion 

We have derived upper limits on the enhancement per unit stellar 
mass rjco of SNIae in globular clusters. The limits are above the 
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Table 5. SNIae for which we were able to exlude a globular cluster origin based on archival HST observations. The time indicates 
the year of the observation used for the exclusion. The VEGA photometric system was used for the apparent magnitudes. CL gives 
the completeness of ruling out globular clusters at the position, see text. 



Supernova Instrument Band Time Observation name Magnitude Galaxy Gal type Distance Ay Type CL 



SN199IT 


ACS HRC 


F814W 


2006 


10607 01 


24.4 


NGC 4527 


Sb 


17.6Mpc- 


0.302' 


P (overturn) 


100 


SN1994aa 


WFPC2 


F791W 


1997 


06419 01 


23.2 


NGC 1320 


Sa 


37.7 Mpc-^"^ 


> 0.2'* 




too 


SN 1 994ae 


WFPC2 


F814W 


2001 


09042 41 


24.8 


NGC 3370 




27.4 Mpc^ 


0.09^' 




100 


SN1994D 


ACS WFC 


F850LP 


2003 


9401 08 


21.8 


NGC 4526 


so 


16.9Mpc2 


0.009' 




50 


SN1998aq 


ACS WFC 


F555W 


2007 


10802 8f 


23.6 


NGC 3982 


s 


22.0Mpc^ 


0.04^' 




50 


SN1998bu 


ACS HRC 


F814W 


2006 


10607 02 


25.3 


NGC 3368 


Sab 


7.2Mpc^ 


0.631' 




100 


SN1999by 


ACS WFC 


F814W 


2004 


9788 13 


25.5 


NGC 2841 


Sb 


14.4Mpc2 


0.030' 


P (sublum) 


100 


SN1999gd 


ACS WFC 


F814W 


2004 


9735 32 


26.1 


NGC 2623 


Sb 


92.4Mpc3 


0.13^' 




75 


SN1999gh 


WFPC2 


F702W 


1999 


06357 14 


24.3 


NGC 2986 


E 


35.3Mpc' 


0.18^' 




75 


SN200()E 


ACS WFC 


F814W 


2003 


9788 d9 


25.1 


NGC 6951 


SBb 


23.0Mpc- 


0.466' 




100 


SN2001el 


ACS HRC 


F435W 


2006 


10883 16 


25.3 


NGC 1448 


Sc 


22.2Mpc' 


0.500' 




50 


SN2002fk 


ACS WFC 


F814W 


2005 


10497 01 


23.5 


NGC 1309 


So 


24.1Mpc2 


0.034' 




75 


SN2005cf 


WFPC2 


F814W 


2007 


10877 18 


24.3 


MCG-0 1-39-03 


SO 


29.9Mpc' 


0.208' 




75 


SN2005el 


WFPC2 


F814W 


2008 


10877 66 


24.3 


NGC 1819 


SBO 


69.5Mpc' 


0.012' 




50 


SN2005hk 


NIC2 


F160W 


2008 


nalp04010 


24.4 


UGC 272 


Sc 


23.0Mpc2 


0.810^ 


P (02cx-like) 


100 


SN2006X 


WFPC2 


F814W 


1996 


06584 02 


23.6 


NGC 4321 


Sc 


13.2Mpc2 


2.496' 




50 


SN2007af 


WFC3 


F160W 


2010 


iblf42010 


22.7 


NGC 5584 


Sc 


22.6Mpc-'' 


0.12-"" 




100 



References: 1-16 see table|3] HBrown et all t20Tot) .^ lTuirv & Fisheil [1981 . ^ IJha et alj (2003) . 



favoured theoretical expectations by almost an order of magni- 
tude, but are well below the observed enhancement of LMXBs, 
iLMXB > 100. We have discussed the effect of the decline of the 
SNIa DTD on the observations. Due to the older age of globular 
clusters compared to the field population in late type galaxies, 
this decreases the value of the parameter rj (which can be seen as 
the average enhancement over the whole population of galaxies), 
which is usually used in the literature (see equation[3]). For this 
reason we find that the sample of SNIae in early-type galaxies is 
more constraining that the full sample, despite containing only 
<l/3 of the supernovae. As this effect of the DTD will be sample- 
dependent we believe that it is more appropriate to use rjco, which 
is the direct measure of the enhancement factor. However, we 
note that it is possible that also rjco is time-dependent, as early 
and late SNIae progenitors are likely to have different evolution- 
ary histories, and since the structure of globular clusters evolve. 

In this pilot study we have only included supernovae which 
had HST images available, and for which we were able to iden- 
tify the position with high accuracy. This sample constitutes only 
a fraction of the total number of SNIae within distances at which 
it is possible to observe globular clusters (~100 Mpc). As no 
globular clusters were detected at the position of the supernovae, 
and as we have managed to probe enough SNIae positions that 
we are starting to constrain theory, it is clearly interesting to ex- 
tend the sample. To reach the values currently favoured by the- 
ory, it will be necessary to expand the sample by a factor of few. 

There are several ways to proceed. As new supernovae are 
discovered continuously, and more galaxies are being observed 
with HST, it is possible to simply wait and let the sample ex- 
pand. However, it is clear that from the current sample that only 
a few SNIae will be added every year, and it will therefore take 
decades before we can reach the goal of a factor of few more 
supernovae. Another option is to use archival data from other 
telescopes. As none of these are as sensitive to globular clusters 
as HST, it will only be possible to probe more nearby SNIae, but 
the advantage is that many more galaxies will be covered by use- 
ful observations. Such a study will be tedious, as data from many 
different telescopes will have to be processed, and it is unlikely 
that a factor of few can be reached. 

A third option is to make a dedicated obsertional program to 
observe the positions of known SNIae. Enough are known that it 
will be possible to reach the goal of increasing the sample by a 



factor of few. The most effective way of reaching it would be to 
make a joint program with a medium-sized ground-based tele- 
scope for the nearest SNIae, and HST for the more distant ones. 
Our experience from this pilot project tells us that it is very dif- 
ficult to obtain useful limits for sources well inside the galaxies 
as the background there is bright and inhomogenous. While the 
number of SNIae will be lower if one excludes sources in the 
inner parts of the galaxies, it will be possible to be complete to 
larger distances, ~ 25 - 30 Mpc with ground-based telescopes 
and 100 Mpc with HST. 

From our results it can clearly be seen that targetting early- 
type galaxies will provide much stronger constraints due to the 
lower difference between the ages (and hence SNIa rate) of the 
field and globular cluster populations, and this effect is being 
enhanced by the higher incidence of globular clusters in these 
galaxies. By specifically targetting SNIae at larger radii, can fur- 
ther reduce the number of observations needed. This is because 
the globular clusters have shallower radial density profiles than 
the stars in the galaxies. For this reason, Fm,cc can be several 
times higher when only considering the outer parts of galaxies. 
It can therefore be possible to probe values of rjco ~few with less 
than 100 positions. 



6. Conclusions 

We have searched for globular clusters at the position of ob- 
served type la supernovae, using archival HST observations and 
literature data. We did not find evidence for globular clusters at 
any of the SNIa positions, and our analysis showed that for 18 
SNIae the observations were sensitive enough that any globular 
cluster should have been detected if it was there. For the posi- 
tions of further 17 SNIae, bright globular clusters would have 
been detected. For the latter sample, we have developed a new 
empirical method to estimate the incompleteness, based on sam- 
ples of globular clusters with bright X-ray sources in nearby 
galaxies. The sample of non-detections allow us to derive up- 
per limits on the fraction of type la supernovae that take place 
in globular clusters. For the full sample the 90% and 99% upper 
limits on the fraction of type la supernovae are 9% and 16%. 
This is higher than the currently favoured estimates of < 1%, 
but within the range allowed by theory. The sample of early- 
type galaxies provides the best limits for the enhancement factor 
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per unit stellar mass for a coeval population rjco ^50 (90% con- 
fidence). We argue that a dedicated survey, combining ground- 
based observations of nearby SNIa positions with HST observa- 
tions of more distant SNIa positions would be able to probe the 
favoured theoretical estimates. 
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Table 6. Upper limits on the fraction of type la supernovae in globular clusters. The table lists the sum of the completeness (CL) 
for each SNIa in the sample. Fj/£_9o% and Fc/i,99% are the upper Umits on the fraction of SNIae in globular clusters derived from the 
sample, < Rqcif > and 77eo,90% are the ranges of assumed values of the fraction of stellar mass in the globular clusters and the SNIa 
rate ratio for the age of the globular clusters compared to the field, respectively. 77co,9()% and 77co,90% are the upper limits on the values 
on the coeval enhancement of the SNIa rate in globular clusters per unit stellar mass, rjco- 



Sample 


Sum 


FuL,90% 


F UL.99% 


Fm,gc 


< Rgc/f > 






Full sample 


27.5 


9% 


16% 


0.3-0.1% 


0.4-0.2 


75-450 


133-800 


Normal SNIae 


21.25 


11% 


20% 


0.3-0.1% 


0.4-0.2 


92-550 


166-1000 


SNIae E+SO 


9.75 


22% 


39% 


0.6-0.3% 


1 


37-73 


65-130 


Normal SNIae E+SO 


7.25 


28% 


49% 


0.6-0.3% 


1 


46-94 


82-163 



